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Abstract 

In this paper we show that the mass-energy density of the Universe can be fully determined in 
terms of the cosmological redshifts, their time drifts and angular-diameter distance (observer area 
distance). Besides providing an indirect measurement of the mass-energy density of the Universe, 
we show how one can use the time-drift of the cosmological redshifts as a replacement for the mass- 
energy density element in the minimally required data set to construct an spherically symmetric 
Lemaitre-Tolman-Bondi (LTB) model for the Universe in observational coordinates. 
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I. INTRODUCTION 



The present astrophysical data sets show that the Universe is almost isotropic around us 
on the largest scales. Using this observed isotropy the Copernican principle leads to a Uni- 
verse modeled by the Friedmann-Lemaitre-Robertson- Walker (FLRW) space-time geometry. 
Relaxing the assumption that the Copernican principle holds, the Universe is then modeled 
by the spherically-symmetric LTB space-time geometry. 

The primary aim of OC program is to strengthen the connections between astronomical 
observations and cosmological theory. We do this by allowing observational data to deter- 
mine the geometry of space-time as much as possible, without relying on a priori assumptions 
more than is necessary or justified. Basically, we want to find out not only how far our ob- 
servable Universe is from being isotropic and spatially homogeneous (that, is describable by 
an FLRW cosmological model) on various length scales, but also to give a dynamic account 
of those deviations (Stoeger et al [lj]). 

Recently, Araujo and Stoeger [2], showed how to construct all spherically symmetric, 
inhomogeneous cosmological (LTB) models in observational coordinates from cosmological 
data on our past light cone, consisting of galaxy redshifts, luminosity distances (observer 
area distances) and mass-energy densities (or, galaxy number counts) as functions of redshift 
and allowing for a nonzero cosmological constant (vacuum energy). In doing so we provided 
a new rigorous demonstration of how such data fully determines the time evolution of all 
the metric components, and a Taylor series algorithm for determining those solutions. This 
enables us to move the solution we obtain from data on the light cone off it in a well- 
defined and straightforward way. It is essential for this to have data giving the maximum of 
the observer area (angular-diameter) distance, Cq{wq, z ma x), and the redshift z max at which 
that occurs. This enables the determination of the vacuum-energy density //a, which would 
otherwise remain undetermined. 

One of the key issues concerning observations to acquire real data is the determination 
of the mass-energy density. For that - if we rely on galaxy number counts - we need an 
accurate determination of the average mass per galaxy as a function of redshift - which 
to be reliable, requires an appropriate model of galaxy evolution which, in turn, requires a 
precise determination of the cold dark matter content in galaxies. 

Here we show that the mass-energy density of the Universe can be fully determined in 



terms of the cosmological redshift, the time drift of the cosmological redshift, hereafter 
redshift- drift, and the observer area distance. Besides its importance in its own right, pro- 
viding an indirect measurement of the mass-energy density of the Universe, thus bypassing 
all the above mentioned problems related to a direct measurement, we can use the redshift- 
drift as a replacement for mass-energy density element in the minimally required data set 
needed to construct an LTB model for the Universe in observational coordinates. Moreover, 
our result is strengthened by noting that the determination of the mass-energy density is 
completely independent of the freedom we have in the observational coordinate formalism 
to specify the time coordinate w along our world line C (see below). 

This relationship of the redshift-drift z to the mass energy-density was first obtained by 



Sandage [3| and McVittie J4], and expressed as 



[l + z)H + H(z) (1) 



lift 



where, Hq is the Hubble parameter now and H(z) is the Hubble parameter at redshi 
z. H{z) is a function of the mass-energy density. Recently, Uzan, Clarkson and Ellis 
working in the 3 + 1 framework, also showed that observation of redshift-drift would allow 
us to test the Copernican principle. The LTB spacetime in the 3 + 1 formalism depends on 
three arbitrary functions of the space-like radial coordinate r, one of which can be fixed by a 
gauge freedom in r. They use this freedom to fix M(r), a quantity related to number counts 
(see below), and proceed to determine the other two using the redshift-drift and observer 
area distance. The results proven here clearly show that in the OC formalism we do not 
need to fix M in that way. Instead, we explicitly determine the mass energy density (number 
counts) as a function of the observer area distance and the redshift-drift. Uzan, Clarkson 
and Ellis 5j also claim that writing the LTB metric in observational coordinates requires the 
solution of the null geodesic equation, which is in general only possible numerically. Given 
the recent results of Araujo and Stoeger mentioned above and those shown in his paper, 
we find that this claim is an overstatement. 

Although we do not yet have the capability to measure redshift drift, the prospects for 
doing so within the next 20 years are good. This should be able to be done using the 
European Extremely Large Telescope (ELT) and its CODEX ultra stable spectrograph (see 
6] and [?j]). Both of these are presently under development. This makes the results given 
here observationally significant. 



In the next section we define observational coordinates, write the general spherically 
symmetric metric using them and present the central conditions for the metric variables. 
Section II I II presents the basic observational parameters we shall be using and several key 
relationships among the metric variables. Section IIVI presents the full set of field equations 
for the spherically symmetric case, with dust and with A ^ 0. In section |V] we present 
our main result showing how the mass-energy density function can be expressed in terms 
of the redshift, the redshift-drift and the observer area distance. In light of this new set 
of data functions, Section I VI I reviews our general integration scheme for all inhomogeneous 
spherically symmetric LTB models. In section IVHI we briefly discuss our conclusions. 



II. THE SPHERICALLY SYMMETRIC METRIC IN OBSERVATIONAL COOR- 
DINATES 

We are using observational coordinates (which were first suggested by Temple 8]). As 
described by Ellis el al [9[ the observational coordinates x l = {w, y, 9, <$} are centered on 
the observer's world line C and defined in the following way: 



(i) w is constant on each past light cone along C, with u a d a w > along C, where u a is the 
4- velocity of matter (u a u a = — 1). In other words, each w = constant specifies a past light 
cone along C. Our past light cone is designated as w = w . 



(ii) y is the null radial coordinate. It measures distance down the null geodesies - with affine 
parameter v - generating each past light cone centered on C. y = on C and dy/du > on 
each null cone - so that y increases as one moves down a past light cone away from C. 



(iii) 9 and are the latitude and longitude of observation, respectively - spherical coordi- 
nates based on a parallelly propagated ortho normal tetrad along C, and defined away from 
C by k a d a 9 = k a d a (f) = 0, where k a is the past-directed wave vector of photons (k a k a = 0). 



There are certain freedoms in the specification of these observational coordinates. In w 
there is the remaining freedom to specify w along our world line C. Once specified there it is 
fixed for all other world lines. There is considerable freedom in the choice of y - there are a 
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large variety of possible choices for this coordinate - the affine parameter, the red-shift, the 
luminosity distance. We normally choose y to be comoving with the fluid, that is u a d a y = 0. 
Once we have made this choice, there is still a little bit of freedom left in y, which we shall 
use below. The remaining freedom in the 9 and (f> coordinates is a rigid rotation at one point 
on C. 

In observational coordinates the Spherically Symmetric metric takes the general form: 

ds 2 = -A(w, yfdw 2 + 2A(w, y)B(w, y)dwdy + C(w, y) 2 dtt 2 , (2) 

where we assume that y is comoving with the fluid, so that the fluid 4- velocity is u a = A~ 1 5% 1 . 

The remaining coordinate freedoms which preserve the observational form of the metric 
are scalings of w and of y: 

~ f dw dy\ 

W^W = w{w), y^y = y{y) I — ^ 7^ — I . (3) 

The first, as we mentioned above, corresponds to a freedom to choose w as any time 
parameter we wish along C, our world line at y = 0. This is usually effected by choosing 
A(w,0). For instance, if the coordinate w is normalized by the condition that it measures 
proper time along C, then A(w, 0) = 1. The second corresponds to the freedom to choose y 
as any null distance parameter on an initial light cone - typically our light cone at w = wq- 
Then that choice is effectively dragged onto other light cones by the fluid flow. Thus we 
here choose y to be co-moving with the fluid 4-velocity, as we have already indicated. We 
shall use this freedom to choose y by setting: 

A(w ,y) = B(w ,y). (4) 

We should carefully note here that setting A(w, y) = B(w, y) off our past light cone w = wo 
is too restrictive. 

In general, these freedoms in w and y imply the metric scalings: 

A^A=^A, B^B = ^B. (5) 
dw dy 

It is important to specify the central conditions for the metric variables A(w, y), B(w, y) 
and C(w,y) in Eq. (j2]) - that is, their proper behavior as they approach y = 0. These are: 

as y -)■ : A(w, y) -)■ A(w, 0) ^ 0, 



B{w,y)-> B(w,0)^0, 

C(w,y)^B(w,0)y = 0, (6) 
C v (w,y)^B(w,0). 

These important conditions insure that C, our world line, is regular - so that all functions 
on it are bounded, and that the spheres (w, y = constant) go smoothly to C as y — > 0. They 
also insure that the null surfaces w = constant are past light cones of observers on C (See 
reference , especially section 3.2, p. 326, and Appendix A for details). 



III. THE BASIC OBSERVATIONAL QUANTITIES 

The basic observable quantities on C are the following: 



(i) Redshift. The redshift z at time wq on C for a co-moving source a null radial distance 
y down our past light cone C~(po), which is centered at our observational position p , is 
given by 

This is just the observed redshift, which is directly determined by source spectra, once 
they are corrected for the Doppler shift due to local motions. It is consistent with and 
complements the first of the central conditions in Eq. OH]). 

We can generalize Eq. ([7|) for the redshift z at an arbitrary time toonC for a comoving 
source a null radial distance y down C~(p) by 

A(w,0) , . 

l + z{ w ,y) = -£-L>. 8 
A(w,y) 

(ii) Redshift-drift (time drift of the cosmological redshift). There are the two fundamental 
four- vectors in the problem, the fluid four- velocity u a and the null vector k a , which points 
down the generators of past light cones. These are given in terms of the metric variables as 

From the normalization condition for the fluid four-velocity, we can immediately see that 
it can be given (in covariant vector form) as the gradient of the proper time r along the 
matter world lines: u a = —r, a . It is also given by fl2]) and OH]) as 

u a = 9abU b = -Aw ta + By >a . (10) 



Comparing these two forms implies 

dr = Adw - Bdy A = r w , B — —r y , (11) 

which shows that the surfaces of simultaneity for the observer are given in observational 
coordinates by Adw = Bdy. The integrability condition of Eq. (fTTj) is simply then 



A' + B = 0. (12) 

where a "dot" indicates d/dw and a "prime" indicates d/dy. This turns out precisely to be 
the momentum conservation equation, which is a key equation in the system and essential 
to finding a solution. 

It follows from Eq. (iTTj) that 

dr ^dw gdy A ^dy 

dz dz dz z dz 

where, z = dz/dw is the redshift-drift. We also note, in passing, that we can write dw/dz = 
dw/dz and dz/dy = dz/dy because w and y are independent coordinates. 

It is clear that on the radial null geodesies r = r(w,y) and z = z(w,y). Now, since on 
our past light cone w = w q - in OC coordinates w labels past light cones - we must have 
that 



dr = k(y)dz 



(14) 



where, k(y) ^ 0. 

Hence, it follows from Eq. (I14p that Jacobian of the transformation {r, z} i— > {w, y} must 
vanish on our past light cone. Since, 



we have that 



J 



d(r, z) 
d(w,y) 



8t dz 

dw dw 

8t dz 

dy dy 



dz 

A{w,y) — 
dy 



B(w,y)z 



(15) 



J(w ,y) = o 



B(w ,y) 



ldz 

z dy 



(16) 



Mwo,y) 

Since z(z) is given from data, solving Eq. ( IT6|) gives z = z(wo,y). That is the same 
information that would be obtained from the null Raychaudhuri [Eq. fl3Tl) below]. It is quite 
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clear that, whereas to solve the later on our past light cone one needs the mass-energy density 
h(wq, z), what we need to solve Eq. f|T6l) is a different piece of information - the redshift-drift 
z(wq,z). Moreover, we shall demonstrate later in this paper that the mass-energy density 
fi mo (z) can be completely determined in terms of the redshift, the redshift-drift and the 
observer area distance C(wq, z) data on our past light cone. 

(iii) Observer Area Distance. The observer area distance, often written as Tq, measured 
at time wq on C for a source at a null radial distance y is simply given by 

r = C(w ,y), (17) 

provided the central condition (jBJ), determining the relation between C(w,y) and B(w,y) 
for small values of y, holds. This quantity is also measurable as the luminosity distance 



because of the reciprocity theorem of Etherington [lO] (see also Ellis 11] ) 



d L = (l + z) 2 C(w ,y). (18) 
(iv) The Maximum of Observer Area Distance. Generally speaking, C(wo,y) reaches a 



maximum C max for a relatively small redshift z max (Hellaby 12]; see also Ellis and Tivon 
[l3| and Araujo and Stoeger [yj]). At C max , of course, we have 



dC{w , z) = dC{w ,y) = Q 
dz dy 

further conditioned by 

d2 °^ z) < 0. (20) 
dz 2 

Furthermore, with the solution of Eq. (Tl6|) . the redshift-drift will give us y — y{z), from 
which we shall be able to find y ma x = y m ax(z max ). These C max and z max data provide 
additional independent information about the cosmology. Without C max and z max we cannot 
constrain the value of A. 

There are a number of other important quantities which we catalogue here for complete- 
ness and for later reference. 

Galaxy Number Counts. The number of galaxies counted by a central observer out to a 
null radial distance y is given by 

N(y) = 4vr f" /i(w , ^m^Biw^ y)C(w , y) 2 dy, (21) 
Jo 



where \x is the mass-energy density and m is the average galaxy mass. Then the total energy 
density can be written as 

dz 1 

K w o,y) = m n(w ,y) = M (z) — — - -, (22) 

ay B(w ,y) 

where n(w ,y) is the number density of sources at (w ,y), and where 

m 1 1 dN 

T tfii rg ~dz' < 23) 
Here dQ is the solid angle over which sources are counted, and J is the completeness of the 
galaxy counts, that is, the fraction of sources in the volume that are counted is J. The effects 
of dark matter in biasing the galactic distribution may be incorporated via m and/or J . 
In particular, strong biasing is needed if the number counts have a fractal behavior on local 



scales (Humphreys et al 15]). In order to effectively use number counts to constrain our 
cosmology, we shall also need an adequate model of galaxy evolution. We shall not discuss 
this important issue in this paper. But, fundamentally, it would give us an expression for 
m = m(z) in Eqs. (1221) and (1231) above. 

The Hubble parameter. The rate of expansion of the dust fluid is 3H = V a u a , so that, 
from the metric (1) we have: 

^(f + # < 24 > 

For the central observer H is precisely the Hubble expansion rate. In the homogeneous 
(FLRW) case, H is constant at each instant of time t. But in the general inhomogeneous 
case, H varies with radial distance from y = on t = to. From our central conditions above 
(3), we find that the central behavior of H is given by 

asy -,0: tfK^^M = tf K 0). (25) 

At any given instant w = wo along y — 0, this expression is just the Hubble constant 
H = H(wo,0) = Aq 1 Bq 1 (B)q as measured by the central observer. In the above we have 
also written A = A(wq, 0) and B = B(wq, 0). 

IV. THE SPHERICALLY SYMMETRIC FIELD EQUATIONS IN OBSERVA- 
TIONAL COORDINATES 



Using the fluid-ray tetrad formulation of the Einstein's equations developed by Maartens 



161 ] and Stoeger et al 



171 ]. one obtains the Spherically Symmetric field equations in obser- 



vational coordinates with A ^ (see Stoeger et al [1] for a detailed derivation). Besides the 
momentum conservation Eq. f fl2|) . they £1X6 £IS follows: 

A set of two very simple fluid-ray tetrad time-derivative equations: 



(B C\ 

u, = -3^ + ^Y (27) 



C\ 6 

where \x m again is the relativistic mass-energy density of the dust, including dark matter, 
and 

1 C C 1 / C \ 2 
^y) = -^ + ACBC + 2{Bc) ' (28) 

is a quantity closely related to ^ mo {y) = fi m ( w o,y) (see Eq l37j) below). 

Equations ff26l) and (j2"T|) can be quickly integrated to give: 

B(w ,y) C 2 {w ,y) 

Hm{w,y) = Vm (y) r 7^7 r; 29 

B(u;,y) C 2 [w,y) 

( fi A \ C 3 (w 0: y) im a 1 C C lfC'\ 2 

u(w, y) = (u (y) + —j - - = -— + — — + - ^— j , (30) 

where co (y) = u(w ,y) and the last equality in (|30|) follows from the definition of u given 
above. In deriving and solving these equations, and those below, we have used the typical 
A equation of state, pa = —fJ-A, where pa and /ia = are the pressure and the energy 
density due to the cosmological constant. Both uq and /io are specified by data on our past 
light cone, as we shall show, /ia will eventually be determined from the measurement of 
C m ax and z max . 

The fluid-ray tetrad radial equations are: 



c 



C \ A B 



(co (y) + f)C 3 (w ,y) 



-\Hmo B(w ,y) C 2 (w ,y) 



c 



+ 



c 



QL — M.QL 

c b c 



(31) 
(32) 
(33) 
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The remaining "independent" time-derivative equations given by the fluid-ray tetrad formu- 
lation are: 



6 CA + A*; (34) 



C CA V 2 

From Eq. ( 152]) we see that there is a naturally defined "potential" (see Stoeger et al [lj]) 
depending only on the radial null coordinate y - since the left-hand-side depends only on y, 
the right-hand-side can only depend on y: 



F(y) = J + (36) 

Thus, from Eq. (1521) itself 

u o(y) = - T^rr 7 / fimoiv) B ( w o,y) C 2 (w ,y) F(y) dy. (37) 

6 2C 6 {w ,y) J 

We now quote a simple but very important observational relationship originally obtained 



by Hellaby 



12] in the 3+1 frameworkjfor a detailed derivation of this result in observational 



coordinates see Araujo and Stoeger [2|] : 

6M max + fi A C 3 max - 3C max = 0. (38) 
where, M(y) is a quantity given by 

M(y) = ±-F mN'(y)F(y)dy = -L f M\y)F(y)dy. (39) 

87T JO 87T JO 

and M max = M{y max ). We note that M(y) = mN{y) is the total mass summed over the 
whole sky by a central observer out to a null radial distance y. 

Eq. (I38p has to be considered a fundamental relation in Observational Cosmology, since 
it enables, from C(wq, z max ) and z max measurements, the determination of the unknown 
constant jjL\. 
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V. MASS DENSITY AS A FUNCTION OF REDSHIFT, REDSHIFT DRIFT AND 
AREA DISTANCE 



Differentiating Eq. flH]) with respect to w gives 



z[w 



A(w,0) { A(w,0) A(w,y) 
V) A{w,y)\A{w,0) A{w,y) 



:i+z) 



A(w,0) _ A{w,y) 
A(w,0) A(w,y), 



Equation ( 153]) can be rewritten as 



(40) 



A_d_(B\ 
Bdw\Aj = 

Substitution of equation (|T2l) into (jHJ) yields 



2u A 2 - -fi m A 2 



(41) 



dw\A ) dw 



l (lnA) 



d_ 

dy 



d_fA' 

dy\A / 



2tu AB + -fi m AB, (42) 



Hence, its general solution is 



A(w,y) = A(w,y)- 



2u(w,y) + -jj m (w,y) 



A(w,y)B(w,y)dy + l(w)\. (43) 



where l(w) is a function of w given by 



l(w) 



A{w,o) 



(44) 



A(w,0) 

It is clear from Eq. (1441) that l(w) depends on the choice of the time coordinate w on our 
world line. 

Substituting Eqs. (EED and $0} into Eq. fl4QD gives 



2u(w,y) + -Hm{w,y) 



A(w,y)B(w,y) dy 



z(w,y) 



(45) 



1 + z(w,y) 7 

It is important to observe at this point that Eq. (j4"5l) is independent of the function l(w), 
thus showing that this result is manifestly independent of any particular choice of the time 
coordinate w. 

Differentiating Eq. ( 1451) with respect to y gives 



2u{w,y) + -fi m (w,y) 



A(w,y)B(w,y) 



d_ 

dy 



z{w,y) 
1 + z(w,y) 



(46) 
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On our past light cone, A(wo,y) = B(wo,y), and we find from Eq. ( H6|) that 



2 d 



1 + zo(y) 



-4w (y), 



(47) 



where, we have denoted A (y) = A(wo,y), and similarly used the same notation for all the 
other quantities. 

Finally, using Eq. ( TIE]) and its solution, we can write Eq. (HTl) in terms of the redshift z 



as 



/Ano \ Z ) 



2io(j a 



Zq(z) 

1 + Zb(*) 



(48) 



Eq. ( T4B]) shows that the mass-energy density fi mo (z) can be completely determined in 
terms of the redshift z, the redshift-drift z(z) and observer area distance C(wq, z) data on 
our past light cone. We observe that the n mo (z) dependency on the observer area distance 
C(wo, z), leads to its dependency on fi\, that must be determined by data giving the maxi- 
mum of the observer area distance, Cq(wo, z ma x), and the redshift z max at which that occurs. 
In the next section we show how this is done. 

Now, from Eqs. ©, flT6j and flU} we find that 



1 



u {z) 



1 + 



z dC (z) 



2C (z) 



, dC (z) 
dz 



(49) 



2Cl(z) y ' Al{z) dz 
Clearly, dCo(z)/dz can be determined from the tq(z) = C(u>o, z) data, through fitting. We 
determine C Q (z) by solving Eq. ([53]) for Co(y) on w = w [see Araiijo and Stoeger [2] for 
details] and then use Eq. ffTB]) and its solution to write the result in terms of z. Taking into 
account the appropriate boundary condition (central condition), Cq(z) is given by: 



Co(z) 



2C (z)Jo z (z) 



Al(z)-zl(~z) d( '" (l) 



dz 



2C (z)dA (z) | dC (z) 
A Q (z) dz 



dz 



A 2 (z)C^(z)j2 A \dz. 



(50) 



Since, z(wo,z) is very small, it follows very clearly from Eqs. ( )48l) . ( 1491) and ( )50l) that a 
useful approximation to fi mo (z) is 



Vm (z) — -4Uq(z). 



(51) 
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VI. THE GENERAL SOLUTION - TIME EVOLUTION OFF OUR LIGHT CONE 



In this section we outline the general integration procedure that is applicable to all inho- 
mogeneous spherically symmetric universe models - that is the only constraint. [See Araujo 
and Stoeger Q for a detailed description of this procedure.] We do not know whether 
the Universe is homogeneous or not. But the data gives us redshifts z, observer area dis- 
tances (angular-diameter distances) r (z), time drift of cosmological redshifts z (z), and 
the angular-distance maximum C max (w ,z) at z max . It is important to specify the latter, 
because, as we have already emphasized, without them, we do not have enough information 
to determine all the parameters of the space-time in the A / case. For instance, although 
we can determine C(w ,z) with good precision (by obtaining luminosity distances d^ and 
employing the reciprocity theorem, equation ffTSj) ) out to relatively high redshifts, at present 
we do not yet have reliable data deep enough to determine C max and z max . But this is just 
becoming possible with precise space-telescope distance measurements for supernovae la. 



In pursuing the general integration with these data, we use the framework and the inter- 
mediate results we have presented in Sections [IV] and |VJ 

First we note that having determined /i mo (y) from data, it is straightforward to obtain 
number counts N(y) from Eq. (jUD - Next, we need the complete determination of the 
"potential" F(y), given by Eq. (136|) . that at this stage is completely determined from data, 
except for its dependency on the unknown So, our next step is to use Eq. (139]) to find 
the mass function M(y). Then we evaluate the mass function M(y) at y ma x and plug the 
result into Eq. ( 1381) which becomes an algebraic equation for fi\. With this determination 
of fj,\, we know Co(y) completely, and can now determine F(y) from Eq. ( 136]) . Furthermore, 
we observe from Eq. (|28|) that the quantity u)o(y) is also completely determined at this 
stage. From here on, we can follow the solution off w = Wq for all w. 

It is shown in Araujo and Stoeger 
following equation for A(w ,y): 



2| that from Eqs. (1121) and (1351) one obtains the 



My) = My)- 



2w (y) + ^ mo {y) 



A 2 (y) dy + C, 



(52) 



where, C\ = A (0)/A Q (0). In the Erratum of reference 2j it is shown that C% and the 
subsequent constants C n that appear in the expressions for higher order derivatives d™A (y) 



14 



(see below) reflect a restricted freedom we have in choosing the time coordinate w on our 
world line. We are free to choose these constants as long as they satisfy the condition of 
being the coefficients in a Taylor expansion of A(w, 0). 

Hence, we have shown that the data on our past light cone determines Aa(y). Since we 
know A (y) from the data, equations (EIJ and f[3"5"j) evaluated on our past light cone become 
algebraic equations for C (y) and B (y), respectively 



h c ( y )My) 
A {y) 

A' (y)A (y) 



<*>o(y) + ^Ad/) 



A 2 (y)C (y) 



A (y) 



2 ^o(y) + q/Wy) 



A 2 (y)B (y) 



(53) 
(54) 



where in the later we have used Eq. (fT2l) . 

Araujo and Stoeger have shown that the next step in this procedure is to obtain and 
equation for A (w , y) through differentiation of Eq. f|35|) with respect to w, and use of Eq. 
(H2]). That leads to: 



( , [v\BS)My) B {y)Al{y) 
A (y) = -A (y) 



{ Al(y) Al(y) 

-^{^[( 2w+ ^)HW- a ° ( ^ 2 ' (55) 

where, C 2 = I (0)/A (0). 

It is important to note that all quantities on the R.H.S. of the above equation are ob- 
tainable either directly from the data or from the algorithmic steps in the procedure we 
are outlining here [see Araiijo and Stoeger |2[ for a detailed description of the algorithm]. 
Therefore, we have shown that we can obtain Ao(y) from the data. The repetition of this 
procedure will give us all time derivatives of A(w,y) B(w,y) and C(w,y) on our past light 
cone, which means that all the metric functions are completely determined by data on our 
past light cone, and calculable as Taylor series. 

We see from the above procedure that each step begins by finding the successive time 
derivatives of the metric function A(w,y) on our past light cone, d™A(wo,y) that is com- 
pletely determined given our choice of the time coordinate w as measuring proper time along 
our world line C. 
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VII. CONCLUSION 



We have shown that the mass-energy density of the Universe can be fully determined 
in terms of the cosmological redshift, the redshift-drift and angular distance (observer area 
distance). This indirect measurement of the mass-density of the universe is instrumental in 
overcoming all the well known and highly discussed problems connected to its determination 
through number counting of sources. That result can also provide a good way of estimating 
the content of cold dark matter in the Universe. 

We have also shown that we can use the redshift-drift as a replacement for the mass- 
energy density element in the minimally required data set to construct an LTB model for 
the Universe in observational coordinates. In fact, this use of the cosmological redshift drift 
as a fundamental piece of data simplifies the integration procedure of the field equations 
in observational coordinates through the use of the much simpler Eq. ([16]) to find z(y) 
instead of the more complex null Raychaudhuri Eq. (l31~j) that in this new scheme, becomes 
a consistency relation. 
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